A model of a pumped continuous atom laser 



Nicholas Robins^, Craig Savage^, and Elena Ostrovskaya^ 

^ Department of Physics, Austrahan National University, ACT 0200, Australia 
^ Optical Sciences Centre, Australian National University, ACT 0200, Australia 



1 Introduction 

At the end of the millennium a continuous wave (cw) pumped atom laser 
remains to be experimentally demonstrated. The prototype atom lasers that 
have been operated have used rapidly pulsed or continuous output coupling, 
but lacked pumping of the trapped atom laser mode , as is needed for a cw 
laser system. Indeed, apart from simple rate equation descriptions, existing 
theoretical models do not incorporate both pumping of the laser mode and 
propagation of the out-coupled beam. 

We have used the macroscopic wave-function approximation, or mean 
field approximation, which leads to the Gross-Pitaevskii (GP) equation, as 
the basis for a complete cw atom laser model. The macroscopic wave-function 
is the order parameter for a Bose-Einstein condensate of atoms. Our model 
system is schematically illustrated in Fig. 0. It includes a pumped reservoir 
of un-condensed atoms irreversibly coupled to the atom laser mode . The 
condensed atoms in this mode, and in the output beam, undergo two-body 
interactions and three body recombination. The atoms are reversibly coupled 
from the atom laser mode to the output beam by a Raman process . We 
solve the resulting system of coupled differential equations numerically in 
one spatial dimension. The macroscopic wave-function of the output atoms 
is analysed to determine quantities such as linewidth. Previous work based 
on systems of coupled equations of the complex Ginzburg-Landau type (GP 
equations with gain and loss) has either neglected important physics, such as 
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Fig. 1. Schematic representation of the model system. The symbols stand for the 
coupling, pumping and loss rates and are defined in the text. 
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the pumping or has concentrated solely on the condensate, ignoring the 
spatio-temporal structure of the output atom beam . 



2 The model 

We consider the trapped atom laser mode, which we refer to as the "laser 
mode" , to be a one-dimensional condensate in a highly anisotropic, cigar- 
shaped harmonic trapping potential. The atoms in the laser mode and in the 
output beam are described by the macroscopic wave-functions Wa(x,t) and 
^b{x,t) respectively. They experience two-body repulsive atom-atom inter- 
actions and three-body recombination. A novel feature of our model is the 
saturated loss due to the three-body recombination process jTj^ . 

To describe pumping we use the phenomenological model of Kneer et.al. 

which mimics the pumping mechanism of a conventional optical laser. 
The number of atoms in the uncondensed pump reservoir is iV„(t), and their 
spatial distribution is not modeled. The reservoir is fed at a rate r, loses 
atoms at the rate juNu, and pumps the laser mode at the rate jpN^Na, 
where Na{t) — \\Pa\^dx is the total number of atoms in the laser mode. 

The laser mode is coupled, by a reversible two-photon Raman transition, 
to an untrapped electronic state which is the output atom laser beam ^bix, t) 
The Raman transition may also give a momentum kick to the out- 
coupled atoms. In an experiment the magnitude of this kick depends on the 
relative geometry of the two Raman laser beams. The output atom beam 
accelerates due to gravity. Inside the trap it overlaps, and hence interacts 
with, the laser mode atoms. This interaction has both a linear contribution, 
due to the Raman coupling, and a nonlinear contribution due to the inter- 
species atom-atom interactions. 

Mathematically, our model is described by the following dimensionless 
equations: 



dt 



r - ^uNu - IpNuNa, 



where Ua and Ub are the intra- and Uab the inter- species two-body interaction 
coefficients, 7^ is the three-body recombination coefficient, and 7^ is the 
Raman coupling coefficient. 
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The model is made dimensionless using the characteristic unit of length 
I = {h/ujmy^'^ and time, r = lu^^, where lu is the trap frequency in the di- 
rection of the weak confinement. We use oj « 125 Hz and m « 10~^^ kg. The 
dimensionally correct two-body interaction coefficients, which we assume 
to be equal, can be written as Ua,b,ab = 47ras/^ where Ug is the s-wave scat- 
tering length for a specific process. There is some arbitrariness in how the 
Ua,b,ab relate to their three dimensional counterparts. The scaling we have 
chosen gives a realistic condensate size Q. Similar reasoning applies to the 
choice of the three-body recombination rate, 7^. We have set it to give feasible 
lifetimes for the laser mode condensate: it is typically between 10~^ and 10~^ 
in our simulations. In accordance with experimentally reasonable values, we 
allow the dimensionless Raman coefficient, = ^2ii72/(wZi), to range up 
to 10^, where ]7i^2 are the Rabi frequencies of the Raman lasers and A is 
the laser detuning Q , and we take the dimensionless Raman momentum kick 
to lie in the range < fc < 100. The dimensionless gravitational coefficient 
G = g{muj~'^fi'^^Y^'^ varies in the range < G < 68, with G = 68 correspond- 
ing to g = 9.8 ms~^. Physically, we may think of the output atomic beam 
propagating in a tilted atom waveguide |9|, with the tilt angle determining 
the value of G. 



3 Rate equations 

Given the complexity of the model equations ([l]), a numerical method is 
the only feasible method of solution. However, some understanding of the 
dynamics of the system can be gained by deriving approximate rate equations 
for the populations of the reservoir, condensate and output fractions JTof . 
We proceed by projecting the condensate and output mean fields onto the 
stationary state modes ^a{x) and 'Pb{x) defined by: 

^a{x, t) = a{t)'^a{x) iTbix, t) = 6(^)^>b(a;)e-^'=^ (2) 

where ^a,6 satisfy the time independent, lossless, and decoupled form of the 
equations (0): 

2 ox^ 2 

2 = -lJ-b<^b + Gx^b + Ub^b\^b\^. 

We are free to choose the magnitude of the chemical potentials //a and /if,. 
In order to simplify the analysis we set the momentum kick, fc, equal to 
zero. By doing so, we neglect the main effect the momentum kick has on the 
dynamics of the system, that is, the change in the spatial overlap between 
the condensate and the out-coupled field. 

The first equation in describes the eigenfunction of a trap with chem- 
ical potential ^a- The equation for the output mode describes a falling wave 
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Fig. 2. Typical population dynamics without three-body recombination, (a) From 
the rate equations (^) showing the growth of the laser mode. We have tested that 
this growth continues in the long time limit (r > 600). Parameters are r = 200, 
7p = 7« = 0.01, 7fl = 0.5, Ua = Ub = 0.02, Uab = 0.01, 5^ = 0.02, Na(t = 0) = 
O.QQl, Nh(t = 0) = 10 and 7r = 0. Values of the spatial integrals are as given in 
the text. The inset shows the Rabi-type oscillations, present in the condensate and 
output fields, (b) From the full GP equations (|l|). Parameters are the same as for 
part (a) and G=68, fc = 5. The GP simulations predict a much smaller number of 
atoms in the output beam. A*';,, than the rate equations, because there is an effective 
loss of atoms due to the boundary absorber at the edge of the numerical grid. 



and can be solved numerically or perturbatively, using a superposition of Airy 
functions in the linear limit We now substitute the ansatz (|]) into the 
original system (|l|), and use the equations (||) to simplify the result. We then 
integrate out the spatial dependence in the equations and separate the time 
dependence of the modal amplitudes as a = no(i)e*^''*^*'' and h = nh{t)e'^^''^*\ 
By noting that the relationship between the modal amplitudes and the pop- 
ulation numbers is given by Na = leSn^ and A^f, = hSnl, we arrive at the 
following system of equations for the real population numbers Na.b, N^, and 
phase difference 9 = Oi, — Oa- 



dt 
dNb 



- {(3Nb + aN^Nl ) - 2-1 R ^hleNaNt sm{9) , 

(4) 

— = 5h + 5nl+ rj^^ cos 61 iVb - Na) 
dt VNaNb 

where S^i = (/i^ - /itf,), a = {2^rh) / , (3 = {2^rh)/{5hY, and 5nl = 

Ub{l - Nb/{5h))l2 + Ua[NJ{5h) - l)h + Uab{{Nah)/{5h) - {Nbl2)/{Sh)). 
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The constant coefBcients Is, where s — 1,2, 3, 4, 5, 6, and S are determined 
by spatial overlap integrals as follows: 



h - J $a^b\^a\^dx, h = J $a^bWb\^dx, (5) 

h = I |<Z'fcpd.T, h = I \'^a\''dx, 5=1 ^a^bdx. 



For realistic parameters, and using the stationary solutions of Eq. (||) we 
numerically calculated these integrals to have the approximate values 6 = 
233, h = 875, h = 1.45, h = 8.39 x 10^ h = 2.79, h = 0.13 and Jg = 41.3. 

From the rate equations it is seen that the populations of the condensate 
and the output mode undergo nonlinear Rabi-type oscillations with the 
effective frequency determined, in part, by the strength of the interactions 
through 5NL{Ua,Uh,Uab), and the Raman coupling 7/j. In the presence of 
pumping and without the three-body loss terms, the condensate fraction ac- 
quires a fast growing component which prevents the system from reaching a 
steady state defined by dNa^b/dt = 0. A typical solution to the system (||) 
is shown in Fig. ^(a). An important feature of this regime is the unbounded 
growth of the laser mode in the absence of three-body recombination, which 
would not be the case for a realistic atom laser. This emphasises the need 
to include three-body recombination in atom laser models, as we do in sub- 
sequent sections. The qualitative predictions of the rate equation model are 
confirmed by direct integration of the system Eq. (|^), with the results shown 
in Fig. |^(b). It is important to note that the values of the spatial integrals, 
Is and 5, play a large part in the time evolution of the rate equations. 
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3. The time evolution of an unpumped condensate. Parameters are 7i{ = 0.5, 
10-^ Ua = Ui = 0.02, Uab = 0.01, G = 12, fc = 5, and Na{t = 0) = 140 
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4 Numerical results 

With this initial understanding of the behavior of the atom laser to guide us, 
we proceed to investigate the numerical solution of our full model equations 
. We use a split-step method to analyse the spatio-temporal behavior 
of the laser mode and the output beam. An absorbing boundary is used at 
the edges of the spatial grid. We have ensured that our results are insensitive 
to the size of the spatial grid and time-step. 

Our model is at its simplest with the pump mechanism turned off, so that 
the out-coupling depletes the atom laser mode. The laser mode population 
changes rapidly at the beginning of the first Rabi oscillation, as can be seen 
in FigJ^, so the approximation of a stationary mode may be false even on 
short time scales One of the results of this paper is that it might be 
more realistic to make this approximation on the assumption of a steady 
state achieved through a balance of pumping and loss due to three body 
recombination. 

We next consider the pumped system. Without the effects of gravity and 
interactions we find the bound state predicted by Hope et al. [^. This is 
produced by the Raman coupling. A fraction of the output field remains 
localised around the laser mode, see Fig.^, and hence the population in the 
laser mode increases indefinitely due to the pumping. Hope et al. found that 
including gravity or two-body interactions could destroy this bound state 
by introducing a repulsive effective potential acting on the output mode. A 
large fraction of the atoms then quickly leave the interaction region. However 
we find that when both gravity and interactions act together they tend to 
recreate the exponentially growing output state. The population dynamics 
of the laser mode and the output is approximately described by Eqs.(|^) in 
this case. The output beam's spatio-temporal structure, in the absence of 
three-body recombination, is shown in Fig.|^. 

Adding three body recombination destroys the bound state as can be 
seen in Fig. || and Fig. 0(a). This occurs because the loss rate increases as the 
square of the atomic density. In addition, the three-body recombination effect 
supresses the atomic density noise in the output field, as well as supressing 
the collective excitations of the condensate, see Fig. 0(b). Remarkably, the 
three body recombination process assists in creating a steady state, both in 
the output and laser modes, as can be seen in Fig. ^{a). 



5 Linewidth narrowing 

One of the characteristic properties of a continuous wave optical laser is the 
dramatic reduction in linewidth that occurs above threshold. The ultimate 
quantum mechanical limit to the linewidth is given by the Schawlow-Townes 
formula . Graham has proposed a related limit for the atom laser due to 

interaction with the thermal reservoir of uncondensed atoms [ p^ . Our model 
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Fig. 4. The bound state of the laser mode (dashed line) and the output beam 
{solid line) at t = 200. The laser mode is multiplied by 0.003. The output beam, 
for X < —2.5 only, is multiplied by a factor of 10. Parameters are r = 200, 7j{ = 0.5, 
7, = 0,Ua = Ub = Uab = 0,G = 0,k = 5, Na{t = 0) = 0.001, Nb{t = 0) = 10. 
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Fig. 5. Atom laser beam density as a function of position and time, with no three- 
body recombination, showing the steady growth of the output. The beam propa- 
gates in the negative x direction from the condensate centered at x = 0. Parameters 
are r = 200, 7h = 0.5, 7, = 0.1, 7p = 0.1, 7, = 0, Ua = Ub = 0.02, Uab = 0.01, 
G = 68, = 5, Na{t = 0) = 0.001, Nb{t = 0) = 10. 
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Fig. 6. Atom laser beam density as a function of position and time, with three- 
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Fig. 7. Steady state operation of the system with three-body recombination, (a) 
Atom number in the output versus atom number in the laser mode. Parameters are 
r = 100, 7fl = 0.5, 7„ = 0.1, 7p = 0.1, 7. = 10-^ Ua = Ut = 0.02, Uat, = 0.01, 
G = 12, = 5, Na{t = 0) = 0.001, Nt(t = 0) = 100. (b) The quasi-steady state 
density profile of the laser mode at t = 150 for pumping of r = 200 (solid) and 
r = 1600 (dotted) Other parameters are 7^ = 0.5, 7u = 0.1, — 0.1, jr ~ 10~^, 
Ua = Ub = 0.02, Uab = 0.01, G = 68, fc = 5, Na(t = 0) = 0.001, Nb(t = 0) = 10. By 
contrast the laser mode at t = 150 for pumping of r = 200 (dashed), with no three 
body recombination, is continuing to grow and change shape. 
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does not incorporate the many-body quantum physics that determines the 
ultimate atom laser linewidth. However, in practice we might expect that the 
linewidth is not determined by fundamental factors, but rather by the dy- 
namics of the particular system; for example, in the optical case by relaxation 
oscillations Our model is suitable for determining the linewidth due to 
this type of dynamical effect. 

To determine the linewidth of the output we have calculated the instanta- 
neous spatial Fourier transforms of particular regions of the output field. The 
results are shown in the inset to Fig.^. The atoms are accelerated by gravity 
according to u = V^o ^ ^.gx and hence, the velocity linewidth narrows as the 
atoms fall away from the trap We are currently investigating the effect 
of the pump rate, r, on the output linewidth of our model, in both the spatial 
and temporal domains. 



Fig. 8. The density of the output beam as a function of position. The inset shows 
the velocity linewidth over three regions of the atom laser output. These regions 
are labeled (a), (b), and (c) and are indicated by the bars on the main curve. 
Their width is about 7.5 dimensionless units, or 68 /im. The number of atoms 
in each region is approximately the dimensionless density times the length; e.g. 
7.5 X 0.4 = 3 for region (b). Parameters are as in Fig.§(a). 



6 Conclusion 

We have presented a model of an atom laser in which we have endeavored 
to include the important physics of an experimental atom laser system. The 
resulting model exhibits a rich and complex dynamics. A novel aspect of our 
model is the inclusion of three-body recombination, which we believe plays an 
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important part in achieving steady state laser operation. We demonstrated 
that the complex interplay of effects in our model can lead to a steady output 
with low atomic density fluctuations. The output of our atom laser is locally 
monochromatic, its linewidth dominated by gravitational acceleration. We 
are in the process of generalising this one dimensional model to three dimen- 
sions. 

We acknowledge continuing discussions with J. Hope, and thank G. Moy 
for his input. 



References 

1. E. W. Hagley, L. Deng, M. Kozuma, J. Wen, K. Helmerson, S. L. Rolston, W. 
D. Phillips: Science 283, 1706 (1999) 

2. I. Bloch, T. W. Hansch, T. Esslinger: Phys. Rev. Lett. 82, 3008 (1999) 

3. G. M. Moy, J. J. Hope, C. M. Savage: Phys. Rev. A 55, 3631 (1997) 

4. J. Schneider and A. Schenzle: App. Phys. B 69, 353 (1999) 

5. M. Edwards, D. A. Griggs, P. L. Holman, C. W. Clark, S. L. Rolston, W. D. 
Phillips: J. Phys. B: At. Mol. Opt. Phys. 32, 2935 (1999) 

6. B. Kneer, T. Wong, K. Vogel, W. P. Shleich, D. F. Walls: Phys. Rev. A 58, 
4841 (1998) 

7. Yu. Kagan, A. E. Muryshev, G. V. Shlyapnikov: Phys. Rev. Lett. 81, 933 (1998) 

8. E.A. Burt, R.W. Ghrist, C.J. Myatt, M.J. Holland, E.A. Cornell, C.E. Wieman: 
Phys. Rev. Lett. 79, 337 (1997) 

9. S. Marksteiner, CM. Savage, P. ZoUer, S.L. Rolston: Phys. Rev. A 50, 2680 
(1994) 

10. N. Robins: Honours Thesis, Australian National University (2000) 

11. J. Schneider and A. Schenzle: App. Phys. B 69, 353 (1999) 

12. R. J. Ballagh, K. Burnett, T. F. Scott: Phys. Rev. Lett. 78, 1607 (1997) 

13. T.R. Taha, M.J. Ablowitz: J. Comp. Phys. 55, 203 (1984) 

14. M. Naraschewski, A. Schenzle, H. Wallis: Phys. Rev A 56. 603 (1997) 



15. J. J. Hope, G. M. Moy M. J. Collett, C. M. Savage: |cond-mat/9907023| , Phys. 
Rev. A 61, 023603 (2000) 

16. A.L. Schawlow , C.H. Townes: Phys. Rev. 112, 1940 (1958) 

17. H.M. Wiseman: Phys. Rev. A 60, 4083 (1999) 

18. R. Graham: Phys. Rev. Lett. 81, 5262 (1998) 

19. A.E. Siegman: Lasers (University Science Books, Mill Valley 1986) 

20. H.M. Wiseman: Phys. Rev. A 56, 2068 (1997) 



